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REAL-TIME  PROCESSING  SYSTEM  FOR  THE 
NRL  REMOTE  SENSING  EXPERIMENT:  PHELPS  BANK,  JULY,  1982 

I.  INTRODUCTION 

This  report  describes  a  software  package  which  was  written  for 
USNS  HAYES  cruise,  July  1982.  The  software  is  a  real-time  processing 
system  which  interfaces  with  the  data  acquisition  software*.  On  the 
HAYES  cruise,  24  channels  of  Meteorological,  Oceanographic,  and  ship 
operating  parameter  data  were  acquired  at  a  rate  of  0.5  Hertz.  The  real¬ 
time  processing  system  facilitated  data  collection  validation,  editing, 
analysis  and  reporting.  The  software  system  operated  successfully  over 
two  eight  day  periods  to  fulfill  the  intended  objectives  described  below: 

1.  Perform  quality  assessment  of  several  data  channels  and 
flag  noisy  data. 

2.  Convert  raw  data  to  physical  quantities  based  on  instrument 
calibrations  and  physical  relationship  between  the  variables. 

3.  Archive  averaged  raw  data  (usually  in  volts  or  frequency), 
calculate  physical  quantities,  and  perform  data  quality  assessments. 

4.  Provide  extreme  execution  flexibility  to  allow  for  accurate 
and  reliable  measurment  of  dynamic  ocean/atmosphere  processes.  The  flex¬ 
ibility  includes  capability  of  on-line  alteration  of  processing  algo¬ 
rithms,  calibration  constants,  and  choice  of  instruments  without  inter¬ 
ruption  of  real-time  acquisition  and  processing. 

5.  Display  in  real-time  the  raw  data  input,  physical  quantities, 
and  measure  of  data  quality. 

6.  Allow  ready  access  to  standardized  file  formats  to  enable 
convenient  and  flexible  off-line  editing  and  analysis. 


*The  data  acquisition  software  was  designed  and  written  by  Dean  Clamons, 
Ocean  Computer  Group  (Code  5003),  Naval  Research  Laboratory. 
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7.  Provide  generalized  processing  techniques  to  permit  future 
use  and  modification  for  tasks  of  a  similar  general  nature. 

The  software  is  disc  resident  on  HP-1000  RTE  VI  system  M2  in  Code 
4310,  Naval  Research  Laboratory,  Washington,  DC  20375. 

After  presenting  a  basic  overview  of  the  data  collection  and  process¬ 
ing  task,  this  report  will  describe  the  real-time  processing  system  from 
three  perspectives: 

1.  Part  III  examines  the  scientific  considerations  of  environ¬ 
mental  instrumentation  and  data  collection. 

2.  Part  IV  promotes  a  user's  understanding  including  instructions 
for  system  implementation. 

3.  Part  V  explains  the  technical  system  software  design  aspects. 
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Figure  2.1.  Fundamental  Processes  of  Data  Collection 


Twenty  four  channels  of  raw  data,  read  by  various  instruments  both 
aboard  ship  and  in  the  surrounding  water,  are  collected  by  the  real-time 
acquisition  program.  Next,  the  data  acquisition  programs  convert  the 
signals  to  digital  form,  archive  the  raw  data  and  pass  the  time  of  day  and 
data  itself  to  the  processing  system.  Data  processing  proceeds  according 
to  guidelines  established  and  maintained  by  the  user.  Every  minute,  data 
averages  are  output  to  disk  and  optionally  to  the  user's  display.  Every  15 
minutes,  data  averages,  results  of  physical  quantity  calculations,  and  data 
quality  assessments  are  output  to  both  disk  and  hard  copy.  The  user  can 
examine  the  output  and,  if  desired,  alter  the  system  execution  guidelines. 


III.  ENVIRONMENTAL  MEASUREMENTS 

The  USNS  HAYES,  July  1982  cruise  was  intended  to  monitor  and  examine 
the  Meteorological  and  Hydrographic  environment  in  shallow  water  during 
radar  overflights  (see  Valenzuela  and  Chen,  1983).  The  software  described 
in  this  report  was  concerned  with  the  time  series  of  meteorological  and 
ship  parameter  data.  A  hydrographic  data  summary  appears  in  Kaiser  (1983). 

A.  Variables  Measured 

The  variables  of  interest  to  the  experiment  were  wind  speed,  wind 
direction,  air  temperature,  atmospheric  dew  point  (or  some  measure  of 
humidity),  water  (bucket)  temperature  and  salinity,  latitude  and  longitude. 
Table  3.1  lists  all  the  input  channels.  Wind  speed  and  direction  measured 
relative  to  the  ship  must  be  converted  to  true  values.  To  do  this  the  ship 
heading(HD),  speed  (SPEED),  and  course  (COURSE)  are  needed.  Further,  the 
wind  sensors  were  situated  on  a  mast  22.5  m  above  the  water  surface.  Ship 
roll  (ROLL)  can  add  a  component  of  motion  to  the  apparent  wind  velocity,  so 
to  correct  for  this,  ship  roll  must  be  measured.  The  USNS  HAYES  is  a  two 
screw  variable-pitch-screw  ship.  The  rotation  rate  (RPM1  and  RPM2)  and 
pitch  of  each  screw  (PI1  and  PI2)  was  logged.  The  ship  also  had  a  pitot 
speed  log  (PLOG)  which  was  logged.  The  latitude  (LAT)  and  longitude  (LONG) 
were  logged  directly  from  a  Northstar  7000  Loran-C  receiver. 

To  measure  these  variables  several  instrument  packages  were 
deployed  on  the  ship.  On  the  starboard  bow  10  m  above  the  water  line,  two 
air  temperature  sensors  (TA1  and  TA2),  an  atmospheric  chilled  mirror  dew 
point  sensor  (TDP1)  and  a  sulfonated  polystyrene  ion  exchange  relative 
humidity  sensor  (TDP2)  were  deployed.  On  the  starboard  of  the  aft  platform 
a  wind  speed  (WS1)  and  wind  direction  (WD1)  sensor  was  installed.  On  the 


TABLE  3.1.  INPUT  TO  PROCESSING  SYSTEM 


Words 

Time 

Description 

0 

Hours 

1 

Seconds 

2 

10’s  Milliseconds 

3 

Years 

Words 

Data 

Channel  # 

Channel 

Description 

0 

0 

TA1 

Air  Temperature  1 

2 

1 

TDP1 

Dew  Point  1 

4 

2 

TA2 

Air  Temperature  2 

6 

3 

TDP2 

Dew  Point  2 

8 

4 

WS1 

Wind  Speed  1,  Starboard 

10 

5 

WD1 

Wind  Direction,  Starboard 

12 

6 

WS2 

Wind  Speed  2,  Port 

14 

7 

WD2 

Wind  Direction  2,  Port 

16 

8 

TS 

Thermosalinograph  Temperature 

18 

9 

CS 

The rmosalino graph  Salinity 

20 

10 

TW 

Water  Temperature 

22 

11 

ROLL 

Ship  Roll 

24 

12 

Spare 

For  future  use 

26 

13 

Spare 

For  future  use 

28 

14 

HD 

Ship  Heading 

30 

15 

RPM1 

RPM,  Starboard  Screw 

32 

16 

PI1 

Pitch,  Starboard  Screw 

34 

17 

RPM2 

RPM,  Port  Screw 

36 

18 

PI2 

Pitch,  Port  Screw 

38 

19 

PLOG 

Pitot  Log 

40 

20 

LAT 

Latitude 

42 

21 

LONG 

Longitude 

44 

22 

SPEED 

Ship  Speed 

46 

23 

COURSE 

Ship  Course 

port  side  of  the  same  platform  WS2  and  WD2  were  installed.  The  water  temp¬ 
erature  (TS)  and  salinity  (CS)  were  measured  in  an  instrument  well,  mid¬ 
ships,  7  m  below  the  water  line.  The  roll  sensor  was  on  the  main 
deck  nearly  midships.  A  towed  water  sensor  (TW)  was  intermittently  de¬ 
ployed  from  the  port  fantail.  Ship  speed  and  course  were  calculated 
and  smoothed  over  about  2  minutes  by  the  Loran-C  receiver  and  then 
logged. 

The  specific  sensors  and  their  characteristics  are 
described  in  Kaiser  and  Munch,  (1983). 

The  sensor  outputs  are  as  follows:  Digital:  TS,  CS ,  PLOG, 

LAT,  LONG,  SPEED,  COURSE.  Frequency:  TW.  Synchro:  HD,  PI1,  PI2.  All 
the  others  are  voltages. 

B.  Processing  Overview 

The  raw  data  was  read  once  every  two  seconds.  This  data  was 
averaged  over  1  minute  to  eliminate  variance  due  to  ship  motion  (typically 
6-12  sec.  period).  Then  various  quantities  were  calculated  from  these 
averages  (correlations  and  variances)  and  these  quantities  were  averaged 
over  15  minutes  to  minimize  the  variability  of  the  lower  atmospheric 
boundary  layer.  The  appropriate  boundary  layer  averaging  interval  is 
actually  determined  by  the  properties  of  the  boundary  layer  existing  during 
the  measurement.  The  appropriate  averaging  period  can  vary  from  100  to 
2000  sec. 

Two  wind  sensors  were  employed.  If  the  wind  was  from  0°  to 
180°  relative  to  the  ship,  WS1  and  WDl  were  entered  into  the  data  stream 
and  if  the  relative  wind  was  from  181°  to  359°  WS2  and  WD2  were  entered 
into  the  processing.  Both  sensor  outputs  were  recorded  in  the  2  second 
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data;  however,  a  user  option  also  allowed  the  selection  of  only  the 
starboard  or  only  the  port  wind  system  data  to  be  entered  into  the  process 
ing  stream. 

The  wind  data  was  corrected  for  ship  roll  since  the  wind 
sensors  were  on  a  high  mast.  Then  the  heading  was  combined  with  the 
wind  direction.  The  wind  velocity  data  was  combined  with  the  ship 
velocity  data  from  the  Loran-C  to  give  true  wind  velocity  at  22.5  m. 

This  wind  velocity  was  then  reduced  to  10  m  elevation  assuming  a  neutral 
atmosphere.  The  TWIND  module  does  the  correction  to  true  wind  at  22.5  m. 

The  dew  point  2  (TDP2)  sensor  actually  outputs  relative 
humidity.  This  is  converted  to  a  vapor  pressure  and  then  to  dew  point 
with  the  use  of  SPOLY.  The  water  temperature  (TW)  sensor  signal  is  a 
frequency  which  was  converted  to  temperature  using  the  manufacturer's 
calibration.  Details  on  the  calibrations,  etc.  can  be  found  in  Kaiser 
and  Munch,  1983. 

C.  Output 

The  original  data  in  digital  format  was  stored  so  that  it  could 
be  rerun  through  the  processing  module  to  reproduce  or  modify  any  calcu¬ 
lations.  However,  the  data  acquisition  programs  did  not  merge  the  Loran-C 
data  with  the  raw  data  and  this  produced  some  downstream  processing  prob¬ 
lems  with  the  2  sec  data  set,  notably  correcting  the  wind  for  ship  motion. 

The  raw  2  sec  data  was  wild  point  edited  and  then  scanned  for 
noise  by  comparing  a  running  60  sec  variance  to  an  average  "background" 
variance.  Any  noisy  records  were  flagged,  and  statistics  on  noise  were 
compiled  for  each  channel.  The  user  selected  the  noise  threshold  for  each 
channel  and  this  threshold  could  be  changed  during  execution. 


The  one  minute  data,  which  represents  atmospheric  conditions 


and  not  ship  motion  contamination,  is  useful  to  examine  variability  on 
short  time  scales. 

The  15  minute  output  is  intended  to  provide  the  final 
characterization  of  the  environment  for  all  the  scientists  participating 
in  the  experiment.  This  output  forms  the  basis  for  the  Meteorological 
data  summary  for  the  July,  1982  experiment  (Kaiser  and  Munch,  1983). 

All  of  the  system  output  is  described  further  in  Chapter  IV  and 


Appendix  A. 


IV.  USER’S  PERSPECTIVE 


A.  Basic  Software  Tasks 

The  real-time  processing  software  consists  of  a  set  of  sub¬ 
routines,  an  interactive  program,  and  a  common  block.  The  "set  of  sub¬ 
routines"  performs  the  actual  data  processing  and  is  controlled  by  the 
calling  module  SAVG§.  The  program  SALTR  allows  execution  characteristics 
of  SAVGS  to  be  altered  interactively  by  enabling  the  user  to  control  exe¬ 
cution  parameters  and  critical  constants.  The  program-to-program  communi¬ 
cation  between  the  SAVGS  modules  and  the  program  SALTR  is  accomplished  via 
SBLK  which  resides  in  a  system  common  area. 

Software  functions  can  be  conceptualized  by  mapping  the  software 
modules  onto  the  Fundamental  Data  Collection  Processes  diagrammed  in 
Figure  2.1,  thus  establishing  the  correspondence  between  necessary  processes 
and  basic  software  tasks. 


Figure  4.1.  Data  Collection  Processes  and  Software  Tasks 


During  real-time  operation,  the  data  acquisition  program  invokes 
a  call  to  SAVGS  every  DT  (2.0  for  this  application)  seconds. 
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During  each  call,  the  data  acquisition  program  passes  SAVGS  the  time  and 
24  channels  of  environmental  data.  SAVGS,  acting  as  the  controlling 
subroutine,  calls  the  remainder  of  the  processing  subroutines  as  needed. 
During  execution,  the  processing  subroutines  are  in  communication  with  the 
common  block  SBLK,  which  contains  guidelines  of  system  execution.  The 
program  SALTR  is  run  by  the  user  when  any  of  the  system  execution  guide¬ 
lines  need  to  be  altered.  The  processing  subroutines  output  the  data 
in  various  forms.  Control  then  returns  to  the  data  acquisition  program. 

B.  Logical  Execution  Flow 

The  next  level  of  description  from  that  shown  in  Figure  4.1 
relates  basic  software  tasks  to  the  system  flow  of  logical  processing 
and  relationship  between  software  elements,  Figure  4.2. 
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In  Figure  4.2,  the  data  averaging  functions,  data  quality  assessment  tasks, 
environmental  calculations  and  output  are  identified.  After  having 
archived  the  time  and  24  channels  of  2  second  raw  data  onto  magnetic  tape, 
the  data  acquisition  programs  pass  the  time  and  raw  data  to  the  processing 
set  of  modules,  specifically  SAVGS.  In  SAVGS,  the  data  is  scanned  for 
extreme  wild  points  and,  if  necessary,  corrected  with  a  magnitude  filter 
The  data  is  corrected  only  to  the  point  that  ensures  continued  program 
execution.  Any  corrected  data  is  flagged  as  it  passes  through  the  software. 
Next,  a  portion  of  the  raw  data  is  converted  to  physical  units,  using  the 
System  Execution  Parameters,  to  allow  for  physical  unit  averaging.  These 
System  Execution  Parameters,  which  reside  in  the  common  block,  SBLK, 
dictate  which  algorithms  are  used  to  calculate  the  physical  units,  the 
channels  to  be  used  in  the  calculations  and  the  instrument  calibrations. 
These  are  discussed  later  in  the  System  Execution  Parameters  portion  of 
this  section. 

One  minute  raw  data  averages  and  physical  unit  averages  are 
formed.  The  time  and  24  channels  of  raw  data  are  written  and  filed  in 
SMIN1.  Also,  using  the  one  minute  raw  data  averages,  data  quality 
assessment  is  implemented  by  running  background  variance  checks.  The 
running  background  variance  algorithm  is  outlined  later  in  this  section. 

If  the  variance  calculations  show  abnormal  results,  a  record  containing 
the  time  of  the  "violation"  and  the  channel  where  the  violation  occurred 
is  written  to  file  SMINSV. 

Fifteen  minute  raw  data  averages  and  physical  unit  averages  are 
formed.  These  15  minute  averages  are  used  to  derive  physical  quantities 
reflecting  environmental  processes.  Data  quality  assessment  is  also 


produced  which  includes  Che  15  minute  total  number  of  background  variances 
and  the  cumulative  number  of  variance  violations  occurring  in  each  data 
channel.  The  15  minute  averages,  calculated  physical  quantities,  data 
quality  assessment  results,  and  the  System  Execution  Parameters  in  effect 
at  the  end  of  the  15  minute  interval  are  output  to  disk  into  file  SMIN5, 
and  hardcopy.  The  hardcopy  15  minute  summary  appears  in  Figure  B.l. 
Descriptions  of  the  disk  files  are  detailed  later  in  this  chapter  in  the 
System  Output  section. 

1.  System  Execution  Parameters 

The  System  Execution  Parameters,  which  reside  in  the  common  block, 
SBLK,  control  execution  of  the  real  time  processing  subroutines  and  are 
set  and  altered  by  use  of  the  interactive  program,  SALTR.  During  normal 
operation  all  parameters  have  to  be  set  before  SAVGS  begins  execution. 
Then,  anytime  after  the  real-time  processing  subroutines  begin  execution, 
SALTR  can  be  run  interactively  to  alter  the  parameters. 

Three  types  of  parameters  can  be  controlled: 

a.  Two  flags  to  indicate  whether  2  second  or  1  minute  averages 
of  raw  data  should  be  monitored  by  dumping  data  to  the  user's  display 
screen  (CRT). 

b.  Values  called  background  variance  thresholds  which  govern  the 
level  of  deviation  allowed  to  be  identified  in  the  19  non-navigation  data 
channels.  If  the  threshold  value  for  a  channel  is  set  low,  there  will 
probably  be  greater  incidence  of  threshold  variance  flags  generated. 

c.  Eight  final  individual  parameters  can  be  viewed  in  two 
different  classes.  The  first  four  control  which  channels  or  procedures  to 


use  during  real-time  processing.  The  second  four  are  actual  values 
involved  in  processing  the  data. 

Choice  or 


Parameter  Default 

1)  TA/TDP  Use  channel  1  or  2?  1  or  2 

2)  Make  roll  correction  to  wind?  (0-No,  1-Yes)  0  or  1 

3)  Make  ship  correction  to  wind?  (0-No,  1-Yes)  0  or  1 

4)  Choose  WS1  to  use  for  wind  calculation?  1 

WS2  to  use  for  wind  calculation?  2 

WS1  and  WS2  *for  wind  calculation?  3 

5)  PA  Atmospheric  pressure  in  millibars  1015.0 

6)  HT  Height  of  air  temperature,  dewpoint,  10.0 

and  relative  humidity  sensors  above 
DWL  (dead  water  line)  in  meters 


7)  HW  Height  of  anemometers  above  DWL  in  meters  22.5 

8)  DT  Time  interval  in  seconds  between  cycles  2.0 

or  scans  of  data  channels 

2.  Background  Variance  Procedure 

A  running  background  variance  check  was  implemented  on  20  of  the  24 
data  channels.  In  the  check,  a  30-point  (60  sec)  running  background 
variance  was  divided  by  the  long  term  variance  for  the  channel  and  the 
ratio  is  then  compared  to  a  threshold.  If  the  threshold  is  exceeded,  the 
channel  is  identified  as  noisy  and  a  record  of  the  violation  is  written  to 
disk  (Refer  to  Figure  A. 2  for  record  format.)  The  thresholds  are 
innitially  set  by  the  user  and  can  be  altered  at  any  time  by  using  the 

*The  wind  direction  dictates  which  sensor  to  use  as  described  in  III.B. 


program  SALTR  to  independently  modify  each  of  the  data  channel  thresholds. 
Total  violations  for  every  15  minutes  are  recorded  on  disk  and  appear  on 
the  15  minute  data  summary  hardcopy  (which  is  essentially  the  same  as 
Figure  B.l).  The  current  background  variance  and  thresholds  for  each 
channel  also  appear  on  the  printout. 

The  algorithm  to  perform  the  variance  check  was  to  be  executed  on 
each  of  the  20  channels  of  data.  Thus,  in  actuality,  during  each 
1  minute  variance  check,  20  variance  ratios  are  calculated  and  compared 
with  20  independent  threshold  values  and  if  necessary,  the  running  back¬ 
ground  variance  must  be  updated.  For  the  first  one  hundred  minutes  the 
background  variance  is  assumed  to  be  a  simple  average  of  the  variance. 

The  algorithm  to  perform  the  variance  check  on  only  one  channel  is  out¬ 
lined  below. 

The  running  variance 


For 


V„  *  1  I  (X.  -  X)2,  X.  -  2  second  data  value. 
b  30  £,!  1  1 

N  <  100,  the  background  variance 


N 


'b  *  t  l  vs 

N  i-i 


where  N  is  the  number  of  minutes  (sets  of  30  data  points). 


For  X  >  100,  Vg  is  calculated  as  follows  (T  is  an  operator  chosen  noise 
threshold) : 


Table  4.1  summarizes  the  output  available  to  the  user. 


Table  4.1.  System  Output 

Disk  Files 

SMIN1  -  contains  time  and  1  minute  averages  of  raw  data 
SMINSV  -  contains  1  minute  data  concerning  data  variance 
SMIN5  -  contains  15  minute  averages  of  raw  data,  physical  data, 
calculated  data  and  variance  data 

Hardcopy 

15  minute  printout,  same  content  as  SMIN5 
Tape  2  second  raw  data 

Optional  user  real-time  display  (controlled  by  setting  System  Execution 
Parameters): 

Time  and  2  sec  raw  data 

Time  and  1  min  averaged  raw  data 

Data  formats  for  disk  files  SMIN1,  SMINSV,  SMIN5  appear  in  Appendix  A. 


D.  Operating  the  Real-Time  Processing  System 

1.  Enable  the  data  acquisition  programs.*  The  major  programs 
are  RDSEN,  RDN7K,  and  STN7K. 

2.  Compile  all  source  codes  for  the  modules  listed  below,  then 
run  the  loader  using  the  following  instructions: 

BGRC 

RE,  ZSENSE 
RE,  %SBLK 
RE,  %SAVGS 
RE,  ZSFILE 
RE,  %SRAW 
RE,  %SAVG1 
RE,  %SCANF 
RE,  %SAVG5 
RE,  %TWIND 
RE,  %SPOLY 
EN 

3.  Compile  and  then  load  the  parameter  alteration  program  SALTR: 
RE,  Z SALTR 

EN 

4.  Run  SALTR  and  fill  in  values  for  all  System  Execution 
Parameters  by  following  the  menu  instructions. 

Terminate  SALTR  at  any  point  after  entering  desired  data. 

♦This  is  done  by  interfacing  with  NRL  Code  5003,  Shipboard  Computer 
Applications  Group 


Start  Data  Processing  Program  ([RU,]  SENSE). 

NOTE:  Make  sure  tape  drives  (LU-8  and  9  usually)  and  the 

printer  are  on  line. 

While  the  system  is  running,  SALTR  can  be  run  anytime  to 
change  the  System  Execution  Parameter  as  necessary. 

To  stop  data  processing  process,  break  the  Data  Acquisition 
Program  (BR,  SENSE). 

NOTE:  All  file  management  is  automatically  performed  as 
long  as  space  remains  on  specified  cartridge  in 
SFILE .  Approximately  50  hours  of  1  minute  data  can 
be  stored  on  one  200  track  cartridge.  About  50  hours 
of  2  sec  data  can  be  stored  on  one  10"-800  bpi  tape. 


V.  SYSTEM  DESIGN 


A.  Module  Calling  Scheme 

A  top-down  modular  design  approach  was  implemented  in  order  to  facili¬ 
tate  overall  task  conceptualization,  scientific/software  functional 
relationships,  and  the  necessity  to  fill  in  software  details  (such  as 
polynomial  instrument  calibration  equations)  as  they  became  available.  A 
hierarchical  tree  structure  having  three  well-defined  levels  was  envisioned 
as  depicted  in  Figure  5.1. 


Figure  5.1.  Hierarchical  Module  Calling  Scheme 


Individual  module  functions  in  the  module  calling  scheme  are 
described  in  the  following  section.  Critical  tasks  such  as  process  syn¬ 
chronization,  dynamic  algorithm  selection  and  mathematical  calculations 
involving  possibly  undefined  results  were  first  tested  in  a  bottom-up 
fashion  before  integration  with  the  top-down  design  structure  to  eliminate 
the  possibility  of  a  real-time  operation  problem.  Finally,  before  field 
operation,  the  entire  real-time  processing  system,  as  shown  in  Figure  5.2, 
was  simulated  using  quality  and  degenerative  test  data. 

System  design  began  in  the  beginning  of  June,  1982  and  was 
completed  in  July  1982  aboard  ship. 


B.  Module  Description 

The  modules  used  in  the  system  are  described  below: 

SAVGS  functions  as  the  main  subroutine,  or  calling  subroutine.  It  also 
implements  a  magnitude  filter  on  19  channels  of  non-navigation 
data.  The  magnitude  filter  clips  data  based  on  preset  maximum 
and  minimum  values  for  each  channel. 

SFI1E  creates  and  opens  output  files  SMIN1,  SMIN5,  and  SMINSV  on  disk. 

SRAW  operates  on  raw  data  using  guidelines  from  the  System  Execution 

Parameters  (see  Section  IV.B.l)  and  acccumulates  2  second  data.  If 
requested,  using  SALTR,  raw  data  can  be  dumped  to  user's  CRT. 

SAVG1  operates  on  1  minute  averaged  data,  writes  1  minute  averages  out  to 
disk,  and  forms  quantities  to  be  used  in  the  15  minute  calculations 
If  requested,  using  SALTR,  one  minute  data  can  be  dumped  to  the 
user's  CRT.  A  running  background  variance  check  is  implemented  in 
SAVG1  to  catch  "noisy"  averaged  data.  One  minute  variances  are  cal 
culated  for  selected  channels  and  if  this  exceeds  the  background 
variance  by  a  predetermined  factor,  then  SCANF  is  called. 

SAVG5  operates  on  15  minute  averaged  data  to  calculate  15  minute  para¬ 
meters.  Results  of  calculations,  data  averaging,  and  background 
variance  data  are  output  to  disk.  A  printout  is  produced  summar¬ 
izing  data  collection  and  analyzed  results  for  the  prior  15  minutes 

TWIND  calculates  the  true  wind  direction.  It  inputs  wind  speed  and 
direction,  ship  speed,  course,  and  heading. 


SPOLY 


calculates  the  polynomial  E(t)  and  T(e).  E(t)  is  the  saturation 
vapor  pressure  of  air  over  water  in  millibars  at  temperature  t. 
T(e)  is  the  temperature  in  degrees  Celsius  at  which  the  atmosphere 
is  saturated  over  water  at  the  vapor  pressure  e  in  millibars. 

SCANF  writes  out  variance  violations  to  disk  and  to  the  user's  CRT. 

SBLK  is  the  common  block.  Extensive  use  of  named  common  is  used  to 
facilitate  modularity.  During  original  use,  SBLK  was  placed  in 
Reverse  Common  as  the  real-time  processing  system  was  implemented 
as  a  large  background  program. 

SALTR  is  the  interactive  program  which  modifies  the  common  block,  SBLK. 

Ultimately  SALTR  controls  algorithm  selection  and  other  parameters 
in  the  SAVGS  group  of  modules. 
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Appendix  A.  Output  Disk  File  Formats 


Table  A.l.  Data  Format  for  File  SMIN1 
(1  minute  averages) 


Words  Channel  #  Channel 

0 

1 

2 

3 


4 

0 

TA1 

6 

1 

TDP1 

8 

2 

TA2 

10 

3 

TDP2 

12 

4 

WS1 

14 

5 

WD1 

16 

6 

WS2 

18 

7 

WD2 

20 

8 

TD 

22 

9 

CS 

24 

10 

TW 

26 

11 

ROLL 

28 

12 

Spare 

30 

13 

Spare 

32 

14 

HD 

34 

15 

RPMl 

36 

16 

PI1 

38 

17 

RPM2 

40 

18 

PI2 

42 

19 

PL0G 

44 

20 

LAT 

46 

21 

LONG 

48 

22 

SPEED 

50 

23 

COURSE 

Description 

Hours 

Seconds 

10' s  milliseconds 
Years 

Air  Temperature  1 
Dew  Point  1 
Air  Temperature  2 
Dew  Point  2 

Wind  Speed  1,  Starboard 
Wind  Direction  1,  Starboard 
Wind  Speed  2,  Port 
Wind  Direction  2,  Port 
Thermosalinograph  Temperature 
Thermosalinograph  Salinity 
Water  Temperature 
Ship  Roll 

For  future  use 
For  future  use 
Ship  Heading 
RPM,  Starboard  Screw 
Pitch,  Starboard  Screw 
RPM,  Port  Screw 
Pitch,  Port  Screw 
Pitot  Log 
Latitude 
Longitude 
Ship  Speed 
Ship  Course 


Table  A. 2.  Data  Format  for  File  SMINSV 
Words  Description 

0  Hours 

1  Seconds 

2  10' s  Milliseconds 

3  Years 

4  Data  Channel  (0-19)  which 

had  a  Variance  Threshold 
Violation 
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Table  A. 3.  Data  Format  for  File  SMIN5  (15  minute  averages) 


Beginning  Word 

Word 

Variable 

Description 

Units 

0 

KD 

Time,  hours  from  beginning  of  year 

Hours 

1 

1(2) 

Time,  seconds  from  beginning  of  day 

Seconds 

2 

1(3) 

10' s  milliseconds  from  beginning  of  second 

10  msec 

3 

1(4) 

Year 

Integer 

4 

LAT 

Latitude 

Deg. 

5 

RLAT 

Latitude 

Min. 

7 

LONG 

Longitude 

Deg. 

8 

RLONG 

Longitude 

Mn. 

10 

LRSPD 

Ship  speed  calculated  by  Loran-C 

Kts 

12 

PLOG 

Ship  speed  from  Pilot  Log 

Kts 

14 

HEAD 

Ship  heading 

Deg. 

16 

CSE 

Ship  course  calculated  by  Loran-C 

Deg. 

18 

TA 

Air  temperature  in  physical  units 

°C 

20 

TDP 

Dew  point  temperature  in  physical  units 

°C 

22 

WS 

Wind  speed  in  physical  units 

m/sec 

24 

WD 

Wind  direction  in  physical  units 

Deg. 

26 

BLANK 

For  future  use 

— 

28 

RLSPD 

Ship  roll  speed 

m/ sec 

30 

RMSRL 

RMS  ship  roll 

Deg. 

32 

TW 

Water  temperature  from  towed  thermistor 

°C 

34 

TS 

Water  temperature  from  thermosalinograph 

°C 

36 

SALIN 

Water  salinity  from  thermosalinograph 

°/oo 

38 

SIGT 

Water  sigma-T  from  TS  and  SALIN 

gm/cm^ 

Table  A. 3.  Data  Format  for  File  SMIN5  (15  min  averages)  (Cont) 


Word 

Variable 

Description 

Units 

40 

TATW 

TA-TW 

°C 

42 

TDPTW 

TDP-TW 

°C 

44 

H 

Sensible  water-to-air  heat  flux 

watts/m^ 

46 

E 

Evaporative  water-to-air  heat  flux 

watts/m^ 

48 

TVATW 

TVA-TVW 

°C 

50 

TVA 

Virtual  air  temperature 

°C 

52 

TVW 

Virtual  air  temperature  when  air  is  saturated 

at  TW 

°C 

54 

RH 

Relative  humidity 

7. 

56 

QWQA 

QW-QA 

gm/kgm 

58 

QW 

Specific  humidity  of  air 

gm/kgm 

60 

QA 

Specific  humidity  of  air  when  saturated  at  TW 

gm/kgm 

62 

TAU 

Momentum  stress  on  water 

dyne/ cm^ 

64 

L 

Monin-Obukhov  length  scale 

m 

66 

ZO 

Sea  surface  roughness  height 

cm 

68 

CD 

Drag  coefficient  for  momentum 

— 

70 

CH 

Drag  coefficient  for  sensible  heat 

— 

72 

CE 

Drag  coefficient  for  latent  heat 

— 

74 

USTAR 

Friction  velocity 

cm/sec 

76 

TSTAR 

TA 

— 

78 

QSTAR 

QA 

— 

80 

PA 

Atmospheric  pressure 

mbar 

82 

HWO 

Height  of  anemometer  above  DWL 

meters 

84 

DT 

Time  interval  between  data  scans 

seconds 

Table  A. 3.  Data  Format  for  File  SMIN5  (15  min  averages)  (Cont) 


Word 

Variable 

Description 

Units 

86 

HT 

Heights  of  air  temperature  &  dew  point 

sensors  above  DHL 

meters 

88 

RBLN2 

Blank 

— 

90 

VB 

Current  background  variances  -  array 

— 

130 

RVAL 

Current  background  variance  thresholds  -  array 

— 

170 

FLG15 

Total  variance  flags  generated  in  last 

15  minutes  -  array 

230 

DAT  15 

15  minute  raw  data  array 

— 

278 

STA 

TA  variance  over  15  minutes 

°C2 

280 

STOP 

TDP  variance  over  15  minutes 

*C2 

282 

SWS 

WS  variance  over  15  minutes 

m2/sec2 

284 

SWD 

WD  variance  over  15  minutes 

deg2 

286 

SPL0G 

PL0G  variance  over  15  minutes 

kts2 

288 

IBLK2 

Blank 

— 

289 

0PTN1 

System  Execution  Parameters  algorithm  selectors 

— 

293 

RPTN2 

System  Execution  Parameters  calibrators 

— 

301 

STS 

TS  variance  over  15  minutes 

°C2 

303 

IBLK5 

Blank 

— 

28 


Appendix  B.  Available  Office  Processing  Programs  to  Use  with  Data  Collected 

SIR  Reads  unedited  1  minute  data  files. 

S5R  Reads  unedited  15-minute  data  files  and  echoprints  data  (Figure  B.l) 
or  writes  a  report  (Figure  B.2) 

35RP  Reads  unedited  15-minute  data  files  and  writes  a  report  which  shows 
results  of  the  15-minute  calculations  (Figure  B.3) 

SNPLT  Plots  2  second  raw  data,  1  minute  raw  data,  and  15  minute  data  and 
calculations.  (Figure  B.4) 

These  processing  programs  are  currently  on  HP-1000  RTE  VI  system  M-2  in 

Code  4310,  Naval  Research  Laboratory. 


HAYES  CRUISE  82-16-09  CC=RM 

192/23S6 

40  43  S3  N  69  17  36  W 

LORSPEED  =  4.4  kts,  PLOG  =  12.  kts  HEAD  =  106  deq  CSE  =  130  deq 

-  CORRECTIONS  TO  10m:  - 

TA  =  14.490  C  TDP  =  14.362  C  WS10  =  S.2?7  MRS  TRUE  UD  =  101  6S 

-  CORRECTIONS  FOR  SHIP  ROLL  - 

ROLL  SPEED  =  .0  m/s  RMS  ROLL  =  6 


TW  =  11  60  C  TS  =  IS  .49  C 

SAL.  =  30.00  SIGT  =  22.02  g/cn**3 


TA-TW  =  2.886  C  TDP-TW  =  2.7S7  C  eW-e  =  -.271E+01  wbar 

II  =  - . 708E+0 1  w/m**2  E  =  - . 346F+02  w/m**2 

TUA-TUW  =  2.91S  C  TUA  =  14  57  C  TMW  =  11. 6S  C 

RH  =  99.1602  X  QW-QA  =  -i.66  q/kg 

GW  =  8.373  g/kq  GA  =10.03  g/kg 

TAU  =  .28  d/cw#*2  L  =  45.4  n 

Zo  =  002873  cm 

Cd  =  92E-03  Ch  =  . 36E-03  Ce  =  .  13E-02 

U*  =  IS . 1 OcM/sec  T*  =  -9E-01  C  Q*  =  .OOE+OOg/kqM 


SIGMAS: 

TA  = 
UD  * 
TS  = 

. 1203E+00 
. 3632E+03 
00  00E+00 

TDP  = 
PLOG  = 
CS  = 

.  1297E  «-00 
.SOiOE+Oi 
. OOOOE+O  0 

WS  =  .  1415E+02 

CHAN. 

IS  FLAG  CUM. 

FLAG  . 

BUAR 

RUAL 

IS  DATA 

0 

0 

0  .  0 

. 1000E-0S 

.  1000E  +  04 

. 14S1E+01 

1 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 1473E+0 1 

2 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 2699E+  0 1 

3 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 4727E *01 

4 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 1379E+  0  0 

S 

0 

0  0 

. 1000E-0S 

.  10  00E+04 

3S04E+00 

6 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 1363E+00 

7 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 3859E  +  0  0 

8 

0 

0  .  0 

. 1000E-0S 

1000E+04 

. 1S49E+02 

9 

0 

0  .  0 

. 1000E-0S 

.  1000E  +  04 

. 3000E+02 

10 

0 

0  0 

. 1000E-0S 

.  1000E  +  04 

. 8221E+04 

11 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

-. 184E+00 

12 

0 

0  .  0 

. iOOOE-OS 

.  1000E+04 

. 1000E+02 

13 

0 

0  .  0 

. 1000E-0S 

.  1000E  +  04 

0000E+00 

14 

0 

0  .  0 

1 000E-0S 

.  1000E+04 

. 1062E+03 

IS 

0 

0  0 

. 1000E-0S 

1000E+04 

12S0E*-01 

16 

0 

0  .  0 

. 1000E-0S 

1  000E+04 

. 3027E+Q2 

17 

0 

0  .  0 

. 1000E-0S 

.  1  000E+04 

2000E-04 

1  8 

0 

0  .  0 

1000E-0S 

.  1  000F+04 

. 6292E+02 

19 

20 

21 

22 

23 

0 

0  .  0 

. 1000E-0S 

.  1000E+04 

. 1266E+02 
. 4081E+02 
-.  69311*02 
.  4422E  ♦•01 

1 301E+03 

OPTIONS 

TA 

TDP=i  ROLL 

COR . WIND 

=1  SHIP  COR 

WIND=i  USi 

. WS2  OR  BOTH 

PA 

=  1013  00 

HT  = 

10.00  HW  = 

22  SO 

DT  =  2.00 

Fig.  B.l.  Sample  o£  the  15-minute  data  and  parameter  summary  provided 


on-line  during  the  experiment. 
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IS  6INUTS  AUF9AGER  DATA 


DTG.Z 

TA 

deqC 

TOP 

deqC 

TS 

deqC 

SAL. 

Pfit 

str,T 

q/CM#*3 

ws 

n/sec 

dsq 

SHIPS 

n/sec 

CS€. 

dt?q 

HDG 

df?q 

l  AT 

dwq  wins 

deq  Mins 

189/0406 

22 

87 

22 

30 

22 

293 

34  071 

23 

448 

4  5 

767 

8 

i 

94 

94 

37 

69 

N 

71 

26 

73 

U 

189/0437 

22 

74 

22 

36 

21 

002 

34  013 

23 

759 

3  2 

266 

6 

9 

92 

95 

37 

51 

N 

71 

18 

27 

U 

139/04S2 

22 

60 

22 

20 

20 

853 

33  898 

23 

712 

3  6 

269 

6 

9 

9.1 

95 

37 

42 

N 

71 

14 

09 

y 

189/0*507 

22 

33 

22 

13 

20 

447 

33.901 

21 

821 

3  9 

260 

6 

9 

91 

96 

37 

39 

”7 1 

9 

84 

j 

189/0522 

22 

01 

21 

94 

20 

246 

33  858 

23 

843 

4  3 

269 

7 

0 

90 

95 

37 

38 

N 

71 

5 

60 

u 

189/0537 

21 

48 

21 

50 

21 

277 

33  879 

23 

583 

4  t 

777 

6 

9 

88 

95 

37 

47 

N 

71 

1 

36 

j 

189/0552 

21 

43 

21 

50 

22 

254 

34  055 

71 

447 

3  7 

278 

6 

4 

80 

95 

37 

73 

N 

70 

57 

25 

y 

189/0607 

21 

86 

21 

92 

24 

998 

35 . 587 

23 

802 

4  0 

255 

6 

5 

84 

96 

37 

1 

14 

N 

n 

53 

89 

y 

189/0622 

22 

28 

22 

34 

24 

147 

35  400 

23 

9\8 

3  6 

270 

6 

6 

36 

96 

37 

1 

35 

N 

70 

49 

S3 

y 

189/3637 

22 

35 

22 

38 

25 

0  37 

35.279 

23 

557 

1  6 

875 

6 

8 

86 

98 

37 

1 

59 

N 

70 

45 

43 

u 

189/0652 

22 

78 

22 

42 

26 

461 

35  839 

23 

534 

4  0 

271 

6 

9 

92 

104 

37 

1 

65 

N 

70 

41 

pp 

y 

189/0707 

23 

46 

22 

29 

26 

R48 

15 . 841 

21 

411 

4  1 

276 

6 

8 

94 

105 

37 

1 

40 

N 

70 

37 

04 

y 

189/0722 

23 

76 

22 

10 

26 

882 
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Fig.  B.2.  Sample  15-minute  average  data  output. 
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Fig.  B.3.  Sample  15-minute  parameter  calculation  output. 
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Appendix  C.  Scientific  Modules  and  Algorithms 


Roll  correction  to  wind: 

HW  =  height  of  wind  sensor  in  m 
DT  =  time  interval  between  data  scans  in  sec 
ROLL  =  roll  sensor  in  degrees  -  current  scan 
ROLLL  *  roll  sensor  in  degree?  -  previous  scan 
WS  M  measured  wind  speed  in  m/sec 
WD  =  measured  wind  direction  in  degrees 
Roll  velocity,  VR: 

VR  -  1.943*HW*(ROLL  -  ROLLL) /(DT*57 .926) 

Across  ship,  WX: 

WMX  =  WS*SIN(WD/57 .296)  -  VR 
WMY  -  WS*COS(WD/57 .926) 

Roll  corrected  wind  speed  WS*  &  direction  WD' 

WS'  =■  [WMX2  +  WMY2]1/2 
WD'  -  arctan  (WMY/WMX) 

Ship  motion  correction  to  wind: 

SPEED  *  ship  speed  in  kts 

COURSE  *  ship  course  in  degrees 

HDG  *  ship  heading  in  degrees 

DIR  =*  wind  direction  in  degrees  (WD') 

SPD  =•  wind  speed  in  knots  (WS') 

TWD  =*  corrected  wind  direction  in  degrees 
TWS  »  corrected  wind  speed  in  knots 
Use  subroutine  TWIND. 
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JOOt  FTN.L 

oaa?  subroutine  tuinocspeed  bourse ,Hor,,orR  ,spd ,tjd , tus> 

0003  D=iHOGtO.TR)  *7  (4159/180  0 

0  004  C=CGURSE*3  14159/130  0 

000S  SUX- SP£fD*C03(C; 

0106  5UY=5PEED*SIN.C> 

0007  twx=5po«cos<d>  -su>. 

0030  TWY=S8D*SIN<D)-Syt' 

0007  TUS*SGRT<  TUX*TUX»TUt*T‘JY> 

0010  TWD=ATAN2(TWY  ,TUX)*lS0  0/3.14159 

0  011  IE<  TWO  LT.0.U>TUD*TUD*.760  0 

0012  RETURN 

0  0 i 7  END 

0014  ENDS 


To  correct  wind  to  10  m  level: 

WS10  *  10  m  corrected  wind  speed 
HW  ■  height  of  wind  sensor  is  a 
WS10  =*  WS/[1  +  .076  ln(HW/10)  ] 

Saturation  vapor  pressure  from  dew  point  or  dew  point  from  saturation 


vapor  pressure: 


IPOLY  =*  1:  RIN  »  Saturation  vapor  pressure  in  mbar 
ROUT  *  Dew  point  in  deg  C 
IPOLY  *  2:  RIN  =  Dew  point  in  deg  C 

ROUT  *  Saturation  vapor  pressure  in  mbar 
The  polynomials  in  SPOLY  were  obtained  by  fitting  data  from  List  (1966) 


0001 
0  0  02 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0  020 
0021 
0  022 
00  23 
0  084 
0025 
0026 
0027 
0  028 
0029 
0030 
0031 
0072 
0033 
0034 


FTN7X,L 

SUBROUTINE  SPOI  Y(DIN,OOUT, IPOLY), 19R2  Cruise  EC 1 3  And  TCe) 


I 


Double  precision  input  To  tuhroutinp 
Double  precision  output  fron  subroutine 
Operational  node  of  subroutine 
Coef  f icien  ts 


Saturation  vapor  pressure  fron  dew  point 


REAL48  DIN 
REAL *8  DOUT 
INTEGER  IPOLY 
REAL48  DO, Di ,02,03,04,05 
C  * 

C  *  E<  t ) 

C  * 

IF  (IPOLY. FQ.l)  THEN 
00*6  1  18348200 
01*0 .4875074800 
D2*l) .  016941.700 
D3*0 . 000119180600 
D4=0  0000061844.700 

DOUT=D0  +  <Di*DIN)  +  (D2*DTN*DIN>  +  <D.7*DIN**3>*<D4*DJN»*4> 
END  IF 
C  * 

C  *  T(e) 

C  * 

IF  < IPOLY. EO  ?>  THEN 
20*-t2. 30400 
01-2.4289300 
02*- 0  062651000 
D3*0  00093679 ’DO 
04*- 0.00 0  0 0689 041  DO 
05*1  936220-8 

OOUT*00  +  (01*OTN)  +  <D2*D tN**2>  +  <D3*DTN«*3>f<  04*DIN**4>+ 
2  <05*OIN**5) 

C  * 

c  * 

c  * 

IF  (OOUT  LT  16  0)  THEN 


Dew  Point  fron  saturation  vapor  pressure 


To  correct  for  hap  fit  of  T(e>  pnlynomal  to  table 
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003S  DOUT  =  0.86<-0  9S400UT 

3036  ELSE 

0037  DOUT=DOUT+0  . 0 1* < DOUT-2?  0)**2-0.31 

0038  END  IF 

0039  C  * 

0040  END  IF 

0041  C  * 

0042  RETURN 

0043  END 


To  calculate  relative  humidity  from  dewpoint: 

E(T)  *  saturation  vapor  pressure  in  mbar  at  temperature  T 

TA  *  air  temperature  in  deg  C 

TDP  *  dew  point  in  deg  C 

PA  *  atmospheric  pressure  in  mbar 

RH  *  relative  humidity  in  percent 

RH  -  100*E(TDP)*(PA  -  E(TA) ) / (E(TA)*(PA  -  E(TDP))) 

To  calculate  dewpoint  from  relative  humidity: 

RH  *  relativity  humidity  in  percent 
R  *  mixing  ratio  over  water  in  gm/gm 
R  -  (RH/100)*.62197*E(TA)/(PA  -  E(TA)) 

E(TDP)  -  R*PA/(. 62197  +  R) 

Use  SPOLY  with  IPOLY  -  1  to  obtain  E(TA)  and  E(TDP). 

To  calculate  sea  water  sigma-T: 


Oj  -  29.42  -  . 270T  -  .0042(T-21)2  +  (.7954  -  .00162T)*(S-34) 

This  is  a  quadratic  fit  to  the  data  in  Stommel,  1965. 

To  calculate  Meteorological  variables: 

CD  *  drag  coefficient 

CD  -  (0.61  +  0.63*WS10)/1000. 

U*  ■  friction  velocity  in  cm/sec 

U*  -  100  /OD*  WS10 

TVA  ■  virtual  air  temperature 

TVA  -  TA(1  +  E(TDP) )  (1  +  .62197  E(TDP) ) 

(PA  -  E ( TDP ) )  (PA  -  E(TDP)) 

TVW  ■  virtual  temperature  of  air  at  water  temperature 

TVW  -  TA  (1  +  E(TW) )  (1  +  .62197  E(TW)) 

(PA  -  E(TW) )  (PA  -  E(TW) ) 

RHO  ■  air  density  in  gm/cm3 

RHO  -  . 34838*PA/ ( 1000* ( TVA  +  273.16)) 

TAU  -  wind  stress  in  dyne/cm2 

TAU  -  RH0*(U*)2 

H  -  sensible  heat  flux  in  watt/m2 

H  -  1046.4*(8.7  +  ,99*(TW  -  TA)*WS10)*RH0 
CH  *  drag  coefficient  for  heat 
CH  -  H/((RH0*WS(TA  -  TW)*1046.4) 

E  ■  evaporative  heat  flux  in  watt/m2 
CE  -  1.32  X  10"3 

E  -  1.515  X  109*RHO*CE*WS*(E(TW)  -  E(TDP))/PA 
Q  ■  specific  humidity  in  g/Kg  at  temperature  T°C 
Q  -  .621*E(T)/PA 
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—  WUUMJV  icuguu  OLdlC  ill  Ui 


L*  • 
p 

L#<  • 


L  =■  267*(TA  +  273.16)*(U*)3*RH0/H 

Zq  *  roughness  height  in  cm 
ZQ  *  (U*)2/79,380 

The  parameterizations  of  CD,  CE,  CH,  E,  H  are  due  to  Smith  (1980). 
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